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1. Introduction. 



In this paper we consider a number of amplitude combinations that are universal, up 
to explicit lattice-dependent factors, in the context of the N ^ limit of the 0{N) model, 
which describes self-avoiding walks and polygons. We correct and extend a table of such 
values given previously by one of us [1], and correct and generalise an exact expression 
for a particular amplitude combination given by the other [2]. We also point out an 
apparent error in another amplitude combination given previously [3]. We make clear 
the role that lattice-dependent factors (such as the number of sites per unit area) play in 
some of these otherwise universal combinations. In particular, we point out for the first 
time a distinction between oriented and unoriented lattices. In order to confirm these 
results, we have generated short series for the radius of gyration of self-avoiding polygons 
on the triangular and honeycomb lattices. We also give arguments for the universality and 
lattice independence of amplitude ratios involving even moments of the area of self-avoiding 
polygons. 

The functions we are considering are: (i) the chain generating function for SAWs, 
C{x) = T,CnX^; (ii) the corresponding polygon generating function, P{x) = T,pnx'"'; (iii) 
the generating function for lattice trails, T{x) = Htnx''^; (iv) the generating function for 
dumb-bell graphs A{x) = T^dnX^; (v) the mean-square end-to-end distance of n-step self- 
avoiding walks (SAWs) {Re)n', (vi) the mean-square radius of gyration of n-step polygons 
{R^)n'i (vii) the mean-square radius of gyration of n-step SAWs {R^)n'i (viii) moments of 
the area of polygons of perimeter n, {aP)n\ and (ix) the mean-square distance of a monomer 
from the origin of n-step SAWs {R^)n- In the above generating functions, Cn,Pnitn and 
dn denote the total number of n-step SAWs, polygons, trails and dumb-bells. For SAWs 
and trails, an origin is chosen on an infinite lattice, and all distinct SAWs and trails are 
enumerated. For dumb-bells and polygons, we adopt the normal convention in which 
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the total number per lattice site is given. As we discuss below, this is different from the 
number of distinct unrooted polygons (up to translations) on certain lattices, such as the 
honeycomb lattice. This distinction is one source of error in Ref. 1. 

In terms of the quantities defined above, we denote the relevant amplitudes as follows: 









Pn 


= BfxV- 






= HfiV 




dn 




'[1 + «(!)] 



(where the growth exponent jj, for trails applies only to lattices of coordination number 
three [4]), and 



{Rein 


= Cn2^[l + o(l)] 


{R\ 


= ^^[1 + 0(1)] 


{aP)n 


= EWn2f'^[l + o(l)] 


{Rg)n 


= Fn2^[l + o(l)] 


{Rm)n 


= G'n2^[l + o(l)] ; 



where 7 = 43/32, « = ^ = g A* = j reciprocal of the critical point. 
Amplitudes -B, D and E^P\ which relate to polygons, will be zero for odd-order terms on 
loose packed lattices, and will be non-zero only for every fourth term on certain lattices 
such as the L-lattice (an oriented square lattice in which every step must be followed by a 
step perpendicular to the preceding step) and the Manhattan lattice. We should remind 
the reader that the above asymptotic forms, and the exact values of the critical exponents, 
are assumptions in the sense that they have not been proven rigorously. Nevertheless they 
all follow from the central assumption that these problems have a continuum limit which 
corresponds to a particular exactly soluble field theory, and it is this correspondence on 
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which shall base our analytic results. The appearance of the same growth exponent n in 
each of the above has been rigorously demonstrated, however. 

In the next section we generalise and correct the first amplitude relation, which gives 
the value of BD. We subsequently discuss the other known relations, and correct the 
relation for the combination BC. We then use these to provide amplitude estimates for 
some cases in which they have not been directly estimated. 
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2. The combination BD. 



Rather than giving the modifications to the argument in Ref. 2 (which in fact contains 
several errors of factors of 2), it is clearer to repeat the whole argument for a general 
lattice. The universality of the combination BD follows from the integral form [5] of 
Zamolodchikov's c-theorem [6], which reads 

c = 37rt'^u-^ J r^{E{r)E{0))cd\ , (2) 

where E{r) is the energy density, which enters the Hamiltonian in the form t J E{r)d\, 
and V is the usual critical exponent governing the divergence of the correlation length. We 
wish to apply this to the 0{N) model whose lattice Hamiltonian is ?i = — a^X^links -^lat(^) 
where r now labels links, and EiQ^^{r) — s{r<) ■s{r>). Here s(r<) and s(r>) are 0{N) spins 
located at the sites r< and r> at the ends of the link at r, ordered in some standard fashion. 
In order to apply Eq. (2) correctly, we must relate the continuum energy density E{r) to 
its lattice counterpart. This is done by equating the continuum and lattice Hamiltonians, 
so that 

t I E{r)d\ {xc -x)Y, ^lat(0 • (3) 

It is convenient to work temporarily on a very large but finite lattice of total area A, and 
to rewrite Eq. (2), using translational invariance, as 

c = 3nt^iy-^A'^ J J^^~ r'f{E{r)E{r'))c d\dV . 
It is then apparent from Eq. (3) that we may simply replace this by its lattice version 

c=M^c-xfv-^A-^Y.^r-r'f{E,,,{r)E,,,{r'))o • (4) 

r,r' 

This should be valid in the scaling limit clS X ^ '^C'l 

when the integral will be dominated 
by values of |r — r'| on the scale of the correlation length ^, which is much larger than the 
lattice spacing, so that the continuum approximation becomes arbitrarily accurate. 
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The next step is to evaluate the right hand side of Eq. (4) in the hmit A?" — > as a sum 
over pairs of mutually self-avoiding walks connecting (r<, r>) and (r^, r>). These pairs of 
walks are then identified with self-avoiding polygons where the links at r and r' are marked. 
The contribution of a given polygon of length n, whose links are at {ri,r2, ■ ■ ■ , rn), to the 
sum in Eq. (4) is then 



where the prime on the sum indicates that adjacent links are to be excluded. This excludes 
contributions where, for example, r> coincides with r^. Such contributions are certainly 
included in the right hand side of Eq. (4), but they are more complicated to evaluate. 
Fortunately they are expected to be negligible in the scaling limit. It is worth checking the 
normalisation in Eq. (5) for small polygons. For example, for n = 4 on a square lattice, 
the sum on the right hand side of Eq. (4) gives 2NJ\fi)X^ where A/5 is the total number of 
links, while Eq. (5) gives ANx^ for each elementary square, of which there are A/j,/2. It 
is straightforward to check this also for larger polygons. Now, if the sum in Eq. (5) were 
unrestricted, it would be equal to 2n^i2^, where is the link- weighted squared radius of 
gyration of the polygon (which differs from the corresponding site-weighted quantity by 
terms of 0(1) as n — > 00.) However, the effect of the restriction is only at the level of 
subleading terms down by one power of n, and is negligible in the scaling limit. If now we 
denote the mean square radius of gyration over all polygons of perimeter n by {R^)n, and 
the total number of such polygons by NsPni where Ms is the total number of sites on our 
large lattice, then, in the limits A?" — > and x — > Xc-, 



This implies that the sum on the right hand side has a singularity of the form {xc — x)~ . 
Since each term in these series is non-negative, this singularity on the positive real axis 
gives the radius of convergence, and, for close-packed lattices, we assume that there are no 




(5) 
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other singularities on \x\ = Xc- Thus this singularity solely dictates the behaviour of the 
terms as n — > oo. For loose-packed and oriented lattices, however, pn is non-zero only if n is 
divisible by an integer a. In this case the generating function is invariant under rotations 
X — > xe^'^^/'^, so that there must exist singularities of equal strengths at a; = Xce^'^^'''/^ 
where r = 0, 1, . . . , cr — 1. Thus, the large n behaviour of pn, for n divisible by cr, is greater 
by a factor a than that expected on the basis of the singularity described by Eq. (6). 



We conclude that 



2^ n C'i^y 



lim npn{R )nX^ = cr^O 



± iv \ - / IV — i\ - — \j r* > 

n— *oo bTT 
where we have introduced the area per site ao = A/Ns- From the known values [7] 
c = 1 — 6/m{m + 1) and i' — (m + l)/4, where N = 2cos(7r/m), for the 0{N) model, the 
second numerical factor is 5/327r^, so that, in terms of the amplitudes defined previously, 
we have 

Bfl=^.a„ . (7) 

For the square lattice, a = 2 and the area per site (measured in units of square lattice 
spacings) is unity. This then gives the result 5/167r^ found previously [2], which agrees 
with results from enumerations given in Ref. 8, as well as later data summarised in the 
final section. 

For a general periodic lattice, we may write the factor ag = A/ Ms as v/k, where v is 
the area of the unit cell, and k counts the the number of sites per unit cell. In the above, 
we have defined pn in the conventional way as the number of n-step polygons per site. If 
instead, as in Ref. 3, it is defined as the number of distinct n-step unrooted polygons (where 
two such polygons are regarded as being equivalent if they are related by translation by 
a lattice vector) then we may associate with each such polygon a unique site. This may 
be done, for example, by ordering the sites in a systematic manner, and associating with 
a given polygon the lowest site according to this ordering. In this way of counting. As 
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should be replaced by the total number of sites which are equivalent to each other by- 
lattice translations. These lattice translations divide the sites into equivalence classes, and 
there is just one representative of each equivalence class corresponding to each unit cell. 
Thus, with this definition of pn, the factor n does not enter the formula for BD. This is 
most easily seen by considering polygons on the honeycomb lattice. The first term in the 
generating function is pQ = ^ if the normalisation is per site, and then k = 2. Alternatively, 
Pq = 1 gives the number of 6-step unrooted polygons. 

This result may easily be generalised to polygons on oriented lattices. In this case each 
link r has a specified orientation, say from r< to r> . Oriented walks on such a lattice are 
described by the N ^ limit of a complex 0{N) spin model, where now the link energy is 
E^^^{r) = s(r>)* • s(r<). Proceeding as before, the sum on the right hand side of Eq. (4) 
is now equal to 

J^E-""' E E'('-t-n)' , 

71 all n— step ^ I 

oriented polygons 

SO that the appropriate lattice-dependent factor is the total number of n-step oriented 
polygons divided by the total area. If we define pn as the number of n-step oriented 
polygons per site, then this factor is equal to Pn/o-Qi as before. However, the complex 
0{N) model is equivalent to a real 0(2N) model, and therefore the appropriate value 
of the central charge on the left hand side of Eq. (6) is c{2N). This leads to an extra 
factor of 2 on diff'erentiating with respect to N in the N ^ limit, so the final result for 
oriented lattices is BD/aaQ = 5/167r^, as compared with Eq. (7) for the unoriented case. 
If we define pn as the number of distinct unrooted oriented polygons, as in Ref. 3, then oq 
should be replaced by v, the unit cell size, in the above formula. 

The results for unoriented and oriented lattices may be unified by supposing that in the 
unoriented case each link may be oriented in either direction, and by defining pn in both 
cases in terms of the number of oriented polygons. The values of pn for unoriented lattices 
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will then be doubled, since each polygon on such a lattice has precisely two orientations. 
The result above will then apply to both types of lattice, and, presumably, to partially 
oriented lattices also. It makes physical sense, because if one looks at polygons with a fixed 
i?^, then the number of oriented polygons on an oriented lattice should be asymptotically 
the same as the number of such polygons on an unoriented, or partially oriented, lattice. 



9 



3. Polygon area amplitudes. 



Next we consider the universahty of amphtudes involving the area of polygons. Our 
approach is as follows: let be an infinitesimal element of the curve formed by the poly- 
gon. Then, by Stokes' theorem, the area within the polygon is proportional io § fx'dl We 
can write this equivalently as §{fx J)dl, where J is a unit current flowing through the 
links of the polygon. Thus, moments of the area distribution are related to integrals over 
correlation functions of J. Such a current will be conserved in the ensemble of polygons. 
In the continuum field theory, it turns out that it is possible to identify the corresponding 
current. Conserved currents in field theory enjoy the special property that their correla- 
tion functions scale according to naive dimensional counting arguments. Moreover, the 
normalization of this particular current is fixed by considering its behavior in the larger 
ensemble including self-avoiding walks with ends. By definition, each end will be a unit 
source or sink for this current. Thus the functional form, and overall normalization, of 
these correlation functions are in principle completely determined in the field theory. 

We now put some fiesh on these remarks. Consider a given oriented polygon, and let 
jl^'t = J^^*) be a unit current fiowing along the links, each of length b and labelled 

by r, in the direction of their orientation. Then the sum 

r 

where ejj is the totally antisymmetric symbol, gives the signed area of the polygon (that 
is, the area is positive for anticlockwise orientation, negative for clockwise orientation.) 
The mean signed area over all polygons with n steps is of course zero, since opposite 
orientations of a given polymer give an equal and opposite contribution. However, the 
mean square area, and, indeed, all the higher even moments, are non-zero and the same 
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as for unoriented polygons. The squared area of a given polygon is 

= E {rir'ijf\r)jf\r') - nr'^jY\r)jf\r')) . 
r,r' 



Using = ?'i</]^*('') = 0, the right hand side may be rewritten as 

Oriented self-avoiding polygons and walks are described by a complex 0{N) spin model 
in the limit A'^ ^ 0. The lattice degrees of freedom are complex 0{N) spins, and the 
Hamiltonian is = —x^j.s{r^)* ■ s(r<) + c.c. Within this model, we may identify 

P^\r) = xb-\r> - f<) (s(r>)* • s{r<) - s{r<r ■ s{r>)) . 

Thus, if (a^)n denotes the mean square area of oriented polygons of n steps, and, as before, 
2J\fsPn is the number of such polygons (the factor of 2 is correct on an unoriented lattice, 
if pn counts the number of unoriented polygons per site), 

'2^fs J2Pn{<^^)nX^ = -H^' E^^ ' r')Hj^^\r) ■ J^^\r')) . (8) 
n r,r' 

— *-| J — * 

In the continuum limit, J is replaced by a conserved current J, which is the Noether 
current for the ^7(1) transformations in the complex 0{N) field theory which commute 
with the 0{N) rotations [9]. The connection between them is simply 



whereby the double sum over r and r' in Eq. (8) becomes a double integral, which, by 
translational invariance, equals 



-\Al? j r^{J{r) ■ J{0))d\ 



(9) 



Now the main point is that the normalisation of J in the continuum limit is fixed by 
the requirement that the unit source for this current is the free end of an oriented SAW [9] . 
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This means that if 0(r) is the continuum field corresponding to s(r) (suppressing 0{N) 
indices), representing the magnetisation density of the 0{N) model, then 

/ J{r') ■ d§r' <^{r) = <^{r) , 

; . (10) 

j J{r') ■ d§r' 0(r)* = -</>(r)* , 

where the integral is over the boundary of a small neighbourhood of r, and d^^r is the 
outward pointing normal. Eq. (10) is supposed to make sense when inserted into correlation 
functions with other operators. It specifies [9] the normalisation of the pole term in the 
short-distance expansion of J with and 0*. 

We see from Eq. (10) that J has unit scaling dimension (since J is a conserved current 
which generates a symmetry of the theory, it has no anomalous dimension.) The fact that 
the normalisation of J is fixed, and therefore does not contain any metric factors, implies 
that its two-point function has the universal form 

( J(r) • J(0)) = rV(r/0 , (11) 

where ^ is the correlation length. Therefore the integral in Eq. (9) is of the form U^'^ ~ 
U^^{1 - x/xc)-^", where U is a universal number, in principle calculable from the under- 
lying continuum field theory, and is the non-universal metric factor in the correlation 
length. 

It then follows from Eq. (8,11) that 

Pn{a\r^aaoU'^^n^''-^x-'' , 

where U' is universal. Now we know that {R^)n is related to the ratio of the second to the 
zeroth moments of the energy-energy correlation function in the 0(N) model, which is of 
the form U"^'^, where again U" is universal. Therefore the amplitude D defined in Eq. (1) 
has the form U"$,q. Also pn ~ Bn~'^^~^x~'^ . Combining these, we find 
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But we argued above that the combination BD/aaQ should be universal and lattice inde- 
pendent. This establishes the universality and lattice independence of the ratio on the left 
hand side, provided the area is measured in square lattice spacings. 

If we now consider the higher even moments of the area, the argument generalises 
straightforwardly. Schematically, 

MsJ^Pni^^^nX^'^hP I ... I ri...rp{J{ri)...J{rp))d\i...d% , 
n J J 

where all indices have been suppressed. The integral on the the right hand side is equal to 

where is universal. Following the argument through, we then discover 

that {aP)n/ {R^)n should be universal, with no lattice dependent factors. 

However, it appears difficult to extend this result to odd moments of the unsigned 
area. This is because these quantities appear not to have any simple expression in terms of 
correlation functions in the continuum field theory. Nevertheless, in Ref. 10 it was argued 
that the corresponding ratio involving the first moment of the area {a})n/ {Fl?)n — E^^^ / D 
should be universal, with no lattice dependent factors. This is to be expected on physical 
grounds [11], since if lattice polygons are regarded as a model for 2-dimensional vesicles, 
the area couples to the pressure p. Since p and (i?^) are macroscopically measurable 
quantities, one might expect a universal relationship between them independent of the 
particular lattice model. Nevertheless, a theoretical demonstration of this is so far lacking. 
In Ref. 12, an argument was made for the finiteness of ratio E^^^/D, but this does not prove 
its universality. Note that [11] the fact that moments of the area satisfy the inequalities 

(a^) < (a2)5 < (a^)5 < ... , 

together with our results for the even moments, does imply that the higher odd moments 
(aP) do scale as expected for p > 3, but it does not imply that their ratios are universal or 
lattice independent. As pointed out [11], it is essential that area be measured in Euclidean 
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units such that the lattice spacing is the same in calculating the radius of gyration and 
the mean area, nd its moments. 
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4. Numerical results. 



Another universal amplitude ratio, first given in Ref. 13 and corrected in Ref. 14 is 

yt,F 2G 1 
yh C C 2 

The individual quotients F/C and G/C are also universal [13]. The best estimates of these 
quantities is given in Ref. 14, as F/C = 0.14029 ± 0.00012 and G/C = 0.43962 ± 0.00033, 
obtained by Monte Carlo calculations. 

The quantity BC/aaQ has been shown in Ref . 3 to be universal. These authors ex- 
pressed their result in terms of BC/va, but they used a definition of pn as the number of 
distinct unrooted polygons. As discussed earlier, these results are completely equivalent. 
It is possible to show directly from the scaling forms of the relevant correlation functions in 
the 0{N) model that C/D is universal. This also follows from Ref. 13 and Ref. 15, where 
it was argued that F/C and F/D, respectively, are universal. Hence, the universality of 
BC/aaQ follows from our result for BD/aaQ. 

The value oi C/D should also be the same for walks on oriented lattices. Since we 
have argued that in this case BD/aaQ gains an additional factor of two, this should carry 
over to BC/aaQ. 

The amplitude for dumb-bells, J, is related to the amplitude for SAWs A by 

J = A{1-2t/ii + t'^/i?)/% , (12) 

as shown in Ref. 16, where r = q — 1 and q is the lattice coordination number. It was 
also shown there that the amplitude for lattice trails, if, on the honeycomb lattice only, is 
related to the SAW amplitude by 

■ (-) 
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Details of the various series and their analysis to estimate critical amplitudes have 
already been discussed in Ref. 1. Since then, extension of the square lattice SAW series 
[17] and the square lattice trail series [18] have allowed more accurate estimates of these 
amplitudes. These are given in Table I below. Corrected amplitudes for the honeycomb 
lattice are also given in Table I. To ensure that the various normalisation factors were 
correct, we generated and analysed the series for the radius of gyratio for triangular lattice 
and honeycomb lattice polygons, and estimated the radius of gyration amplitude directly. 
These new series are given in Table II, and the amplitude estimates in Table I. The quantity 
(triangular) was given in Ref. 1 in units of unit triangles. It is more correctly given in 
Table I in units of area, assuming unit lattice spacing, as also used in the estimate of other 
metrically dependent amplitudes. Direct estimates of the amplitudes A, B and C for the L 
and Manhattan lattices are also given. These have been obtained from the series in Ref. 19 
and Ref. 20. The polygon counts for Manhattan lattice polygons quoted in Ref. 19 should 
be divided by 2 in order that the correct normalisation per-site be retained. 

The square lattice amplitudes are generally the most accurate. We see that 
327r2SL>/5(Tao = 1.0000 for the square lattice, 1.020 for the honeycomb lattice and 0.997 
for the triangular lattice. This is perfect agreement, given the accuracy of the amplitude es- 
timates, and allows -D(L, Man.) to be estimated. Similarly, the invariant BC/aao = 0.2168 
(square), 0.2167 (triangular) and this value is used to predict C(honeycomb). The ratio 
F/C = 0.1403(square), 0.1402 (triangular), in good agreement with the Monte Carlo es- 
timate cited above. This then permits F(honeycomb, L, Man.) to be estimated. We see 
that this quantity is a factor 2 greater for the oriented L and Manhattan lattices, as pre- 
dicted. This factor was not found by Privman and Redner [3] , who found that (with their 
convention for pn as the number of distinct unrooted oriented polygons) the same result 
for BC/va as for unoriented lattices, within their numerical errors. It seems likely that 
these authors miscounted by a factor of 2. For example, the number of distinct oriented 
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4-step polygons on the L lattice is 2, not 1. 

Similarly, G/C = 0.4397(square), and 0.4402 (triangular), again in agreement with the 
precise Monte Carlo estimate. The predicted values of G(honeycomb, L, Man) are given 
in Table I. The ratio E/D (where E is E^^^ in our previous notation) was found to be 
2.515(square) and 2.529(triangular), where we believe the square lattice estimate to be 
more precise, as the series from which E was estimated is far longer for the square lattice. 
This value then permits (honeycomb, L, Man.) to be estimated. The amplitudes J 
follow from the amplitudes for A and Eq. (12). The amplitude (honeycomb) follows from 
A(honeycomb) and Eq. (13). All amplitudes quoted are expected to have an associated 
error confined to the last decimal place quoted. Apart from some minor gaps for the 
triangular, Manhattan and L lattice amplitude. Table I now gives a complete and corrected 
tabulation of critical amplitudes. 

Acknowledgments. We would like to thank Prof. Michael E. Fisher for valuable comments 
and correspondence, and Mr. Andrew Conway for programming assistance in calculating 
the radii of gyration. This work was supported in part by the Isaac Newton Institute, the 
UK Science and Engineering Research Council, the US National Science Foundation Grant 
PHY 91-16964, and by the Australian Research Council. 



17 



Tables 



Table I. Estimates of critical amplitudes A ... J defined in the text, for the honeycomb, 
square, triangular, L and Manhattan lattices. Quantities in parentheses are estimates 
from the amplitude relations discussed in the text. The remainder are estimates from the 
corresponding series expansions. 



Amplitude 


Honeycomb 


Square 


Triangular 


L 


Manhattan 


A 


1 

± . ±^0 


1 1 771 

-L . -L / / -L 


J- . -LOU 


-L .yjo 


n SQ 

\J .OiJ 


T3 
J_> 


u.uooo 


U.OOZO 




9 47 


9 1^ 


c 


(0.889) 


0.771 


0.711 


0.67 


0.73 


D 


0660 


05631 


0.0518 


(0.049) 


(0 053) 


E 


(0.166) 


0.1416 


0.131 


(0.12) 


(0.13) 


F 


(0.125) 


0.1082 


0.0997 


(0.095) 


(0.10) 


G 


(0.389) 


0.339 


0.313 


(0.30) 


(0.32) 


H 


(1.341) 


1.272 








J 


(0.000972) 


(0.002768) 


(0.006205) 






1^ 


1.8477591 


2.6381585 


4.1507951 


1.5657 


1.733 


a 


2 


2 


1 


4 


4 


ao 


3^3/4 


1 


V3/2 


1 


1 


3277"^ BD/5aao 


1.020 


1.0000 


0.997 






BC/aao 




0.2168 


0.2167 


0.41 


0.46 
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Table II. Number of polygons and the sum of the squares of the radii of gyration of n-step 
polygons on the triangular and honeycomb lattices. 



Triangular 




Honeycomb 


n 


Pn 


Ain?pn{B?)n 




n 


'iPn 


Sn'^Pn{B?)n 


3 


2 


24 




6 


1 


144 


4 


3 


96 




8 








5 


6 


408 




10 


3 


2460 


6 


15 


1872 




12 


2 


3168 


7 


42 


8688 




14 


12 


32052 


8 


123 


39912 




16 


18 


77976 


9 


380 


183264 




18 


65 


420444 


10 


1212 


834744 




20 


138 


1310088 


11 


3966 


3779064 




22 


432 


5655204 


12 


13265 


17013936 




24 


1074 


19291968 


13 


45144 


76186320 




26 


3231 


76066992 


14 


155955 


339566616 




28 


8718 


268063080 


15 


545690 


1507025568 




30 


25999 


1011675420 


16 


1930635 


6662739096 










17 


6897210 


29355291552 










18 


24852576 


128932421592 
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